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Abstract. We present an algorithm for computing discriminants and prime ideal decom- 
position in number fields. The algorithm is a refinement of a p-adic factorization method 
based on Newton polygons of higher order. The running-time and memory requirements 
of the algorithm appear to be very good: for a given prime number p, it computes the 
p-valuation of the discriminant and the factorization of p in a number field of degree 1000 
in a few seconds, in a personal computer. 



1. Introduction 

The factorization of prime numbers in number fields is a classical problem, whose resolution 
lays at the foundation of algebraic number theory. Although it is completely understood from 
the theoretical point of view, the rising of computational number theory in the last decades 
has renewed the interest on the problem from a practical perspective. In his comprehensive 
book [CohOOi p. 214], H. Cohen refers to this problem as one of the main computational tasks 
in algebraic number theory. 

The most common insight in the known solutions of the problem is based on the solution 
on a more general problem: the determination of a (local) integral basis. There is a number 
of highly efficient methods for this problem, due to H. Zassenhaus and M. Pohst |PZ89| . D. 
Ford and P. Letard |FL94j . and D. Ford, S. Paufi and X. Roblot |FPR02| . 

The theory of higher order Newton polygons developed in [Mon99j and revised in [HNj . 
HN standing for "higher Newton" , has revealed itself as a powerful tool for the analysis of the 
decomposition of a prime p in a number field. Higher Newton polygons are a p-adic tool, and 
their computation involves no extension of the ground field, but only extensions of the residue 
field; thus, they constitute an excellent tool for a computational treatment of the problem. In 
this paper we explain how the theoretical results of [HNj apply to yield an algorithm, due to 
the second author [Mon99[ Ch.3], to factor a prime number p in a number field K, in terms 
of a generating equation f{x). The algorithm computes the p-valuation of the index of f{x) 
as well; in particular, it determines the discriminant of the number field, once one is able to 
factorize the discriminant disc(/) of the defining equation. 

In many applications, the computation of an integral basis is very useful because it helps to 
carry out other tasks in the number field. However, if one is interested only in the discriminant 
or in the factorization of a prime, our direct method has the advantage of being more efficient 
and it makes possible to carry out these tasks in number fields of much higher degree. In fact, 
the running-time and memory requirements of the algorithm appear to be very good. Even in 
some bad cases, chosen to test the limit of its capabilities, it computes the factorization of p 
in a number field of degree 1000 and p-index 200000 in a few seconds, in a personal computer. 
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If we add the computation of generators of the prime ideals, the running-time may increase 
in a significant way, because this routine imphes an extended gcd computation. 

The outUne of the paper is as follows. In section 2 we present the basic algorithm that is 
obtained by a direct application of the ideas of [HNJ . In section 3 we introduce an optimization 
based on a lowering of the order in which the computations take place, and we prove a strong 
optimization result (Theorem 13. ip . We refer to this optimization process as refinement, and 
it results in a dramatic lowering of the complexity. In section 4 we show how to compute 
generators of the prime ideals lying above p in terms of the output of the algorithm. In section 
5 we describe an implementation, and in section 6 we present the results of some numerical 
tests. We construct some "worse possible" polynomials, that should be specially difhcult with 
respect to the structure of the algorithm; this means that they have a huge index, and this 
index is sufficiently "hidden" to force the algorithm to work in a high order. The record is a 
polynomial of degree 6912 and 2-index 77673504, for which the factorization of 2 is obtained 
in 787 seconds. The algorithm, moreover, is highly parallelizable, so that it can raise the 
bounds of computations on number fields to huge degrees. 

The local nature of all the computations involved in the algorithm justifies its high efficiency 
compared to the classical insight explained above. Anyway, after this algorithm, one can go 
the other way round and apply it as a previous step in the determination of an integral 
basis. Numerical experimentation suggests that this new approach provides a significant 
improvement in the solution of this problem. 

2. Computation of discriminants and prime ideal decomposition in number 

FIELDS 

We fix a number field K — Q{0), generated by a monic irreducible polynomial f{x) £ Z[x], 
such that f{6) — 0. We denote by Zk the ring of integers of K. We fix also a prime number 
p e Z. The p-adic valuation is denoted simply by v in order to avoid confusion with p-adic 
valuations Vr of higher order. If F is a finite field and <p(y), "0(2/) S '^[v]: we write ~ -0 to 
indicate that the two polynomials coincide up to multiplication by a nonzero constant in F. 

In this section we present the basic algorithm that computes the p- value of the discriminant 
of K and the prime ideal decomposition of p , that is obtained by a direct application of 
the ideas of [iM] . 

2.1. Types ands their representatives. The basic tool for the algorithm is the concept of 
type and its representative, which we recall here with some detail. All results of this section 
are taken from [HNl §2]. 

Definition 2.1. A type of order zero is a monic irreducible polynomial in ¥p[y]. Let r > 1 
be a natural number. A type of order r is a sequence of data: 

t {(l)i{x); Xi,<j>2{x); . . . ; Xr-i,(j>r{x); Xr,ipr{y)), 

where (f)i{x), . . . , 4'r{x) G Z[a:] are monic polynomials, Ai, . . . , G are negative rational 
numbers, and ipriv) £ '^piu] monic polynomial, that satisfy the following properties: 

(1) 4>i{x) is irreducible modulo p. Let ipo{y) € Vp[y] be the polynomial obtained by reduc- 
tion of <j)i{y) modulo p. We define Fi :— Fp[y]/('(/'o(y))- 

(2) For all 1 < i < r, the Newton polygon of i-th order, Ni{(j)i^i), is one-sided, with 
positive length and slope Aj. 

(3) For all 1 < i < r, the residual polynomial of i-th order with respect to Xi, Ri{(j)i^i){y) , 
is an irreducible polynomial in Vi[y\. Let ^i{y) G ^i[y\ be the monic polynomial 
determined by Ri{(j)i+i){y) ~ tpi{y). We define F^+i ¥i[y]/{')pi{y)). 

(4) ^priy) G ^r[y] is a monic irreducible polynomial, tpriy) 7^ y- We define F^+i := 

Vr[y]/{A{y))- 
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Every type carries implicitly a certain amount of extra data, whose notation we fix now. 
For all f < i < r: 

• hi,ei are a pair of positive coprime integers such that = —hi/ci, 

• ii, &Z are fixed integers such that iihi — i'^a — I, 

• fi^ degV'i(y), and /o = degV'oly) = deg0i(a;), 

• TTii — deg0i(a;), and m^+i = nirerfr- Note that m^+i = eifiUii, 

• z^^y (mod V»(y)) S F*+i, zq = y (mod Vo(2/)) £ F^. Thus, F,+i = F,(z,). 

Also, for all 1 < i < r+1, the type carries certain p-adic discrete valuations Vi : Qp(a;)* ^ Z 
[HN| Def.2.5], and semigroup homomorphisms, 

Lo,: Iplx] \ {0} ^ Z>o, P(2;) ^ ord^,_, (A,_i (P)), 

where Ro{P){y) G Fp[j/] is the reduction modulo p of P{y)/p^^-^\ These objects play an 
essential role in what follows, because Ldi{P) measures the length of the principal part, N~{P), 
of the Newton polygon of i-th order of P{x) [HN , Lem.2.17]. The principal part of a 
polygon N is the polygon determined by the sides of negative slope of N. 

To avoid confusion, in case of working simultaneously with different types, we add a su- 
perscript with the type to every component or datum: '/'*(a;), A*, e\, etc. 

Definition 2.2. We say that Ai, 0i+i(a;) (and their implicit data) are the i-th level oft. 

By truncation we can easily obtain types of lower order. We denote 

t; := (0i(x); Ai, (/)2(a;); • • • ; A^-i, 0i(a;); Ai, V'i(y)), < i < r, 

ii := ((/)i(a;); Ai,02(a;); ■ ■ • ; Ai_i, ^^(a;); A^, 0i+i(a;)), < i < r. 

Clearly, tj is a type of order i. The extension ti is not a type, strictly speaking. 

To our polynomial f{x) G Z[a;] we can attach the set to(/) of all types of order zero that 
divide f{x) modulo p. By Hensel's lemma, each t e to(/) determines a monic p-adic factor 
ft{x) G Zp[x] of f{x), and 

/(^)= n /t(^)- 

teto(/) 

The types of order r play an analogous role and they provide similar factorizations in higher 
order. Let us recall some concepts and results in this respect. 

Definition 2.3. Let t, t' be types of order r, and let P{x) G I'plx] be a monic polynomial. 

• We say that P{x) has type t if degP = mr+iUJr+i{P) > 0, or equivalently 

(1) P{x) = 0i(x)"° (mod p), for some positive integer oq, and 

(2) For all 1 < i < r, the Newton polygon Ni(P) is one-sided, of slope Xi, and 
Ri{P){y) ^ Il^iiy]; for some positive integer a^. 

• We say that P{x) is divisible by t, or that t divides P{x), if uj^_^i{P) > 0. Formally, 
we can think o/ a;*_|_]^(P) as the exponent with which t divides P{x). 

If t divides P{x), we denote by Pt{x) the monic factor of P{x) of type t and greatest 
degree. It has degPt = ™r+i^r+i(^); ^'^d uj^_^i{Pt) = u}^j^i{P). 

• We say that t and t' are P-equivalent, if both divide P{x), and Pt{x) ~ Pt'{x). 

• We say that a set T of types faithfully represents P{x), if P{x) is divisible by all types 
in T, and P{x) — IlteT^tl^)- 

In [HNl §2.3] it is described a constructive method to enlarge a type of order r into different 
types of order r + 1 . 

Theorem 2.4. Let t be a type of order r. We can effectively construct a monic polynomial 
ipr+iix) G Z[a;] of type t such that Wr+i(0r+i) — 1- This polynomial has minimal degree 
deg0r+i = "T-r+i among all polynomials of type t. 
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We call such a polynomial (j3r+i{x) a representative of the type t. We denote by t := 
{(j>i{x); • • • ; Ar, (l)r+i{x)), the extension of t by (j)r+i{x); this object is prepared to be enlarged 
to a type of order r + 1, (t; A^+i, V'r+ily)), simply by taking any negative rational number 
Xr+i G and any irreducible monic polynomial ipr+i{y) € ¥r+i[y]. 

The representative of a type plays the analogous role in order r to that played by an 
irreducible polynomial modulo p in order one. 

2.2. Types versus prime ideals. The Basic algorithm. Recall that we have fixed a 
monic irreducible polynomial f{x) G 1i[x]. 

Definition 2.5. A type t of order r is said to be /-complete if ujl_^^(f) — 1. 

Theorem 2.6 ( |HN1 Cor. 3. 8]). Lett be an f -complete type of order r. Then the p-adic factor 
/t(x) is irreducible in l^plx]. Moreover, if L/Qp is the extension generated by ftix), we have 

e{L/Qp) = ei • • • e,, f{L/Qp) = /o/i •••/.. 

Thus, an /-complete type singles out a unique prime ideal p dividing pZx, whose ramifi- 
cation index and residual degree can be read in the data of t: 

e{p/p) = ei • • • e^, /(p/p) = /o/i • • • /r- 

The p-adic factorization process of |HN| consists essentially in the construction of a set T 
of /-complete types, that faithfully represents /(x). Thus, it can be interpreted as a Basic 
algorithm, to determine the prime ideal decomposition of pZj<-. The types are built iteratively 
by means of Theorem 12.41 and the theory of Newton polygons of higher order. We start with 
the set To(/) := to(/), that faithfully represents f{x). We extend the non- /-complete types 
of this set to types of order one, in order to construct a set Ti(/) that, again, faithfully 
represents f{x), etc. At each order r, the extension process is carried out by a main loop that 
performs the following operations. 

Main loop of the Basic algorithm. At the input of a non- /-complete type t of order 
r — 1, for which LOr{f) > 0, and a representative (firix): 

1) Compute the Newton polygon of r-th order, Nr{f) = Si + ■ ■ ■ + St, with respect to t 
and 4>r{x)- 

2) For every side Sj of negative slope A^j < 0, compute the residual polynomial of r-th 
order, Rr,j{f){y) G II^r[2/], with respect to t, (l)r{x) and A^.j. 

3) Factorize this polynomial in [y] : 

4) For every factor ?/'r,fc(y), compute a representative of the type t^''' :— (t; \rj,tpr,k{y))- 
For those factors 4'r,k{y) with exponent Ok — 1, the type t^''' is complete. For the remaining 
types we continue the iterative process. 

Thus, each non-complete type of order r — 1 has sprouted several types of order r, which 
are called branches of the input type t. We have a factorization in Zp[x]: 

ft{x) = J|/t3.fc(2;), 

with deg/t3,fc — Cr.j fr,k'mr- Also, (ojr+i)*^ (/) > 0, for all j, k, and 
(2.1) io'rif) = ^e,,,/,,fcK+if' '(/). 

Hence, the invariant Ltj*{f) is an upper bound for the number of irreducible factors of ft{x), 
and it is a kind of measure of the distance that is left to complete the analysis of the type 
t and its branches (or equivalently, to decompose each fl'^(x) into a product of irreducible 
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factors). Also, (12. ip shows that, except for the case in which there is only one branch with 
Cr = /r = 1, the branches are always closer to be /-complete than t. 

We denote by tr{f) the set of types of order r obtained by aplying the Main loop to all 
non- /-complete types of tr_i(/). We denote by ti{fy°™P^ the subset of the /-complete types 
of ti(/), and we define 

T,(/):=t,(/)U [ U t,(/r-P>| . 

Proposition 2.7 ( HN, §3]). Tr{f) faithfully represents f{x). 

To show that the Basic algorithm deserves this name, we have to prove that, after a finite 
number of enlargements, all types of tr(/) will be complete. To this purpose we introduce 
another variable to measure how far is a type from being complete, that works even in the 
unibranch case with = /r = 1. This control variable is defined in terms of higher indices. 

2.3. Indices of higher order. The results of this section are extracted from |HN[ §4]. De- 
note 

ind(/) -vHIk- m)). 
and recall the well-known relationship, v(disc(/)) = 2ind(/) + v{disc{K)) , between ind(/), 
the discriminant of f{x) and the discriminant of K. 

Definition 2.8. Let N = Si + ■ ■ ■ -\- St be a principal polygon, with finite sides ordered 
by increasing slopes Xi < ■■■ < Xt < 0. Denote by Ei = i{Si), Hi ~ H{Si), dt = d(Si) 
the respective length, height and degree of each side |HNl §1.1]. We define the index of the 
polygon TV to be the nonnegative integer 

* 1 

ind(iV) ^ E,H,. 

i=l l<i<j<t 

This number is equal to the number of points with integral coordinates that lie below or 
on the polygon, strictly above the horizontal line that passes through the last point of N and 
strictly beyond the vertical axis. Hence, ind(iV) = if and only if N has a unique side with 
height H — 1, or length E — 1. 

Definition 2.9. Let t be a type of order r — 1, and let 4ir{x) be a representative of t. We 
define its /-index to be the nonnegative integer 

indt(/) :=indt.^^(/) := /o • • • /._i ind(7V-(/)), 

the Newton polygon of r-th order taken with respect to t and (j)r{x). 
We say that t is /-maximal ift divides f(x) and indt(/) = 0. 
For any natural number r > 1, we define indr(/) := X^tGt i(/) i^dt(/)- 

Since the Newton polygon N~{f) depends on the choice of (prix), the value of indt(/), and 
the fact of being /-maximal, depends on this choice too. 

Proposition 2.10 ([HNI Lem.4.16]). 

a) Ift is f -complete, then it is f -maximal. 

b) Ift is f -maximal, then either t is f -complete, or the output of the Main loop applied 
to t is a unique branch of order r + 1, which is f -complete. 

Thus, the fact that all types of tr{f) are complete is essentially equivalent to the fact that 
they are all maximal. The proof that this will occur after a finite number of iterations is 
provided by the Theorem of the index. 
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Theorem 2.11 (Theorem of the mdex pHl Thm.4.18]). For all r>l, 

(2.2) ind(/)>indi(/) + ... + ind,.(/), 

and equality holds if and only if all types oftr{f) are f -maximal. 

This theorem shows that after a finite number of iterations aU types of tr{f) will be /- 
maximal, because the sum of the right-hand side is bounded by the absolute constant ind(/). 
By (b) of Proposition 12.101 either Tr(/), or Tr+i(/), will contain only /-complete types. 
By Theorem 12.61 and Proposition 12.71 our family of complete types determines the complete 
factorization of pTix into a product of prime ideals. At this final stage we have necessarily 
an equality in (|2.2p . so that we get a computation of ind(/) as a by-product. 

Remark 2.12. If at the end of step 1 of the Main loop of the Basic algorithm, we accumulate 
to a global variable the value indt(/), the final output of this global variable is ind(/). In 
particular, indt(/) is an absolute measure of the distance covered by each iteration of the 
Main loop, towards the end of the algorithm. 

Summing up, we have proved the main theorem of the paper. 

Theorem 2.13. Given a number field K, a generating equation f{x) G 1j[x], and a prime 
number p, we can construct a set T of f -complete types, that faithfully represents f{x). The 
types of T are in 1-1 correspondence with the prime ideals of K lying above p, and the 
ramification index and residual degree of each ideal can be read from data of the corresponding 
type. Along the construction of T, the algorithm computes the p-valuation of the index of 
f{x) as well. 

The Theorem of the index and Proposition 12 . 101 show that the number of iterations of the 
Main loop is bounded by ind(/). Actually, in practice, the number of iterations is much lower, 
because in each step, indt(/) is usually much bigger than one, due to the abundance of the 
number of points of integer coordinates below an average Newton polygon with a fixed length 
LJ*(/), and the fact that these points are counted with weight /o • • • fr-i- 

In the next section we introduce a crucial optimization. A refinement process will control 
at each iteration wether it is strictly necessary to build a type of higher order, or it is possible 
to keep working in the same order, to avoid an increase of the recursivity in the computations. 
For instance, the polynomial f{x) = {x — 2)^ -I- 2^*^ would require the construction of types of 
level w fc in a strict application of the Basic algorithm, while it can be completely analyzed 
with a refined type of order 1. 

3. Optimal representatives of types 

3.1. Detection of optimal representatives. The construction of types dividing a given 
polynomial is not canonical: in the construction of the representatives (j)r(x) one makes some 
choices, mainly related to lifting certain polynomials over finite fields to polynomials over Z. 
A natural question arises: are there some choices better than other ones? 

Consider the following trivial example: let p = 2, f(x) = x'^ — Ax -\- 12, and K = Q{d) — 
Q(\/^), with Zk — The polynomial f{y) has only one irreducible factor, ^po{y) = y, 

modulo 2; thus, the type of order zero t = 4'oiy) gives no information about the factorization 
of 2 7iK- The more natural lifting of V'o to Z[a;] is 4>i{x) = x, and the corresponding Newton 
polygon and residual polynomial determine a unique extension of t to a type of order one, 
{x;—l,y + 1), which is still not complete, so that we must construct a type of (at least) 
order 2 to determine the factorization of 2Zx- If we choose instead, 4>i{x) = x — 2, we find 
f{x) — {x — 2)^ -{-2^, and the unique extension of t to a type of order one, (^i(a;); — 3/2, ?/ + 1), 
is complete. Thus, it is clear that this second choice of (f>i{x) is better. 
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While it seems very difficult to predict a priori whether a choice of (j)r{x) is better than 
another, it is possible a posteriori to know if our choice was optimal and, if this is not the 
case, to improve its quality. 

Theorem 3.1. Let t'^ e tf.-i(/) be a type of order r — 1, which is not f -complete, and let 
4>r{x) he a representative of t'^ . Let t — (t'^; A^, '0r(?/)) S '^rif) be one of the branches of t'^ , 
and let ft{x) G Zp[a;] be the factor of f{x) determined by t. Let (j)'j.{x) G Z[X] be another 
representative oft^. Lfcrfr > 1, then, 

a) The Newton polygon N^(ft), with respect to t*^ and (f)'^{x), is one-sided with slope 
'^'r ^ '^r, o,rid it has the same end point than Nr{ft)- 

b) The residual polynomial i?J,(/t )(?/), with respect to t", 't'ri^) '^"'^ K' ^'^^ '^''^^V ''^^ 
irreducible factor in Vr[y]; that is, R'^{ft){y) ~ ^^(y)"'-, for some monic irreducible 
polynomial '0r(y)- 

c) Let (y)' = ; A^-i, 0r(a;)) he the extension o/ t° determined by the choice 
of(t)'^{x), and let t' ((tO)'; A;, V-^Cy))- Lf > Xr, then = /; = 1. // A; = A^, 
then = Cr, f^ = fr, and (^l+iif) = ^t+iif)- 

Therefore, if e^/r > 1, the representative (j)r{x) is optimal for this branch t of order r. The 
absolute measures indto^^/ (/), indto are not the right invariant to compare because they 

incorporate the influence of other branches. If we center our attention on the branch t, there 
are two situations in which the choice of (/>5,(a;) would lead us to be closer to the end of the 
analysis of this branch: 

(1) t is replaced by several branches of order r, 

(2) t is replaced by t', with uj^^^if) < 

Items a), b) show that any choice of (j)r{x) leads to replacing the branch t of t°, by a single 
branch t'. Also, if AJ, = A^, we have uJ^_^_l{f) — Ljl^i{f); thus, replacing the type t by t' 
makes no difference at all in this case. 

However, if |AJ,| < |Ar|, we get a definitely worse approximation to the final solution, 
because = {f)/{erfr), by (|2.ip . Thus, the type t is much nearer to be complete 

than t'. Also, if > 1, indt(/) will be probably bigger than indt'(/), because each point 
of integer coordinates below N~_^j^{f) will contribute with a higher weight, /o ■ ■ ■ fr, to the 
/-index. 

Note that a choice of the representative (prix) of t° can be optimal for some branches t 
and non-optimal for other branches. We shall see later that the condition e^/r > 1 is also 
necessary for the optimality of (prix) with respect to the branch t. 

For the proof of Theorem l3.11 we need an auxiliary result. Fix a type t" of order r—1 and di- 
viding f{x). For any n — {uq, . . . , rir-i) G N*", denote ^{n){x) = . . . (pr-iix)^'—^ . 
Let 6* G Qp be a root of /to(a;), and L ^ Qp{9). In [HNl (27)], an embedding Fl, is 
defined by 

(3.1) Le : ¥r ^ F^, zq 9, zi ^ • ■ • , z^-i ^ jr-i{0), 

for certain rational functions 7i(a;) G Z(x) such that v{'ji{9)) = [HNi Def.2.13,Cor.3.2]. 

Lemma 3.2. Let t'^, 9, L be as above. Let M{x) G T^x] be a polynomial of degree less 
than TUr- Suppose that n G N*" satisfies v{M{9)) — v{^{n){9)). Then, the nonzero element 
M{9)/^{n){9) G F^ belongs to teiWr), and the element Lg^{M(9)/<^{n)(9)) G F* is indepen- 
dent of the choice of 9. 

Proof. Let J ■— {i = {jo, ■ ■ ■ , jV-i) G N'' | < ji < e^/i, for < i < r}, where we take eo = 1 
by convention. Since degM < mr, this polynomial can be written in a unique way as 

Mix)^ ajx^''$(0,ji,...,>_i)(x). 
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for certain integers aj. By |HN[ Lcni.4.21], wc have 

v{ai) > Si v{M{d)) - «(<i>(0, ji, . . . 

for all j e J. Let Jq = {i ^ J \ v{aj) = Sj}. Denote 6j = cij/p^', and j' = ((Jj, ji, ■ • ■ ,jr-i)- 
We can write M(x) as 

M{x) = bix^'>m'){x) + N{x), 

jGJo 

where N{x) e Z[x] satisfies v{N{9)) > v{M{9)). Now, 

*(n)(a;) ^ ■> $(n)(x) 

Jfc Jo 

By hypothesis, - n)(6l)) = 0. Since - n)) = [Ml Prop.2.15], we have 

^'r($(j'-n)(x)) = 0, by |HNl Prop. 2. 9]. By [HN) Lem.2.16], there exists a sequence zi, . . . , V-i 
of integers, that depend only on j' and n, such that 

m' ^n){x)=ji{xy' ■■■jr^xY'-K 

Hence, the element of F^, 

belongs to Lg{¥r). Since all the ingredients aj, Sj, ii,. .. ,ir-i etc. depend only on t*^, the 
element L~\M{9)/U{d)) G ¥; is independent of 6. □ 



Proof of Theorem \3.1[ Let 9 £ Qp he now a root of ft{x), and L — Qp{9). Let us prove first 
that v{(j)r{0)) > v{4i'r{9)). In fact, let us show that v{4)r{9)) < v{(f)'^{9)) implies e,- = fr = 1- 
Let us write (p'^ix) — <j)r{x) + M{x), for certain polynomial M{x) S Z[a:], of degree less than 
mr. Ifv{(f>ri9)) < w(^;(6')), then v{M{9)) = vl(f>ri9)) = + |A^|)/ei---e^-i, by the 

Theorem of the polygon |HN| Thm.3.1]. Since degM < rrir we have LUr+i{M) = |HN| 
Lem.2.2], and [HNl Prop.2.9] shows that Vr{M) = ei • • • er-iv{M{9)) = Vr{4>r) + |Ar|; hence 
Xr is an integer, and Cr = 1. 

We use now some other rational functions introduced in [HNi Def.2.13], and the identity 

Vr{(l)r) — Gr-lfr-lVr{4>r-l) |HN1 Thm.2.11]: 

_ ^r{x) _ (j)r{x) 

~ T:r{xy'^ ~ ^.(x)'''-7r,_i(x)''-(^-)/«-i ■ 

Denote n(a;) = TTr{xf^iTr-i{xY-'^'l'-^l''--K Since v{-ir{9)) = 0, we have v{n{9)) = v{(j)r{9)) = 
v{M{9)). By |HN[ (17)], we can write W.{x) = ^{n){x), for some n G N'' that depends only 
on t". Now, if we reduce modulo ttil the identity 

n(^)) -^"^^^'^ n{9))' 



Lemma 13.21 shows that 7r(^) = —M{9)/Il{9)) belongs to ig{¥r)- Since 7r(^) is a root of 
ieiAiy)) [HNl Prop.3.5], we get = 1. 

We prove now a) of the theorem. If we show that v{4)'^{9)) takes the same value for all the 
roots 9 of /t(x), then, by the Theorem of the polygon, N^{ft) will be one-sided with slope 

= Vr{(t>'r) - ei • • • er^iv{(t>'^{9)) = Vr{(t>r) - Ci ■ ■ ■ er-iv[(f)'^{9)) 

> Vr(4>r) — ei • • • er~-iv{<pr{9)) = \r- 

Now, if v{(l)'^{9)) = v{(j)r{9)) for all 9, then the value v{(j)'^{9)) is constant, because the 
value v{(j)r{9)) is constant. Note that v{M{9)) — Vr{M)/ei ■ ■ ■ e^-i is independent of 9. 
If there is one 6*0 with v{<p'J9o)) < 'y(0r(6'o)), then v{M{9o)) = w(0;(6'o)) < viM^o))- 
Hence, v{M{9)) < v{(t)r{9)) for all 9, because both expressions are independent of 9. Thus, 
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v{(j)'j.{9)) — v{M{9)) is constant too. Finally, the polygons Nr{ft), N^{ft), have both end 
point iLufif),vf{f)). 

We prove items b), c) simultaneously. Suppose first that > A,.; then, the Theorem 
of the polygon shows that v{(j)'j,{6)) < v{cj)r{9)). Arguing as above, this implies = 1 and 
7^(6') e te(Fr), with r/ := ig^{j'r{(^)) ^ independent of 9. By the Theorem of the residual 
polynomial |HN[ Thm.3.7], if 9 runs on all the roots of ft{x), then 7^(0) runs on all the roots 
of the irreducible factors of i?^(/t)(y). Hence, R'^{ft){y) ^ {y — '7)"'', and = 1. 

Suppose now AJ, — A^, so that v{(t)r{9)) — v[(j)'^{9)), by the Theorem of the polygon. 
We distinguish two cases. If v{M{9y) — v{(pr{9)), arguing as above we get = 1 and 
lr{9) = 7r(^) + ie('7), for some 77 € F* which is independent of 9. In this case, i?J,(/t)(y) ~ 
Rr{ft){y + 77) is a nonzero constant times the power of an irreducible polynomial of degree 
/; = fr. If v{M{9)) > v{(f >ri9)), t hen ^ 9^'- = (ti^^Y'- + N{x), where v{N{9)) > viM^Y")- 
Arguing as above, we get jri9) = ^^{9), and this implies R'^{ft){y) ~ Rr{ft){y) and /,', = fr- 
This implies = t^l+iU) too, because erfrU^t+iU) = e'J^uj^'^^if) by dHH). □ 

3.2. The process of refinement. What can be said when Cr = fr = 1? In this case, we 
enlarge the type t*^ to an order r type t — (t°;— — 77), and we find a representative 
4>r+i{x) of t, of degree ttj^+i — rrir. Let us emphasize a crucial observation. 

Remark 3.3. The polynomial (t>'r{x) ■— (f>r+i{x) can be taken too as a representative of t*^. 

In fact, (t)r{x) has type t°, and iOr{<j)'r) — deg(j)'r/mr — 1. We shall show that (j)r{x) is 
always a better representative of t" than 4>r{x); thus, in this case 4'r{x) is never optimal. 

The comparison betwen these two representatives is done by means of the following affine 
transformation: 

n-.M.'^ — H{x,y) ^ ix,y - hrx). 
Note that the vertical lines of the plane are invariant under this transformation, and Ti. acts 
as a translation on them. Also, Ti. preserves points of integer coordinates. If S' is a side of 
negative slope, of length £, slope A and degree d, then 7i(5') is a side of length £, slope X — hr 
and degree d. 

Definition 3.4. Let h be a positive integer, t a type of order r — 1, 4>r{x) a representative of 
t, and P{x) G Z[x] a nonzero polynomial that is not divisible by 4>rix). 

If N is a principal polygon, we denote by the part of N formed by the sides of slope 
less than —h. We define 

(3.2) md'liP) := ind^^^^jP) := /o • • • fr-i {\nd{Ni\P)) - hi{e ~ 1)^ , 

where I — £{Nr{P)), and the Newton polygon is taken with respect to t and (prix). 

This number indt(P) is equal to /o • • ■ fr-i times the number of points of integer coor- 
dinates in the region of the plane determined by the points that lie below (or on) Nr{P), 
strictly above the line L^h of slope —h that passes through the last point of the polygon, 
and strictly beyond the vertical axis. The term hl{£ — l)/2 takes care of the points of integer 
coordinates in the triangle determined by L-h, the vertical axis and the horizontal line that 
passes through the last point of iV/'(P). 

Let us introduce some notation. Let iVr+i(— ), denote the Newton polygon with respect to t 
and (j)r+i{x). Let N^{—) denote the Newton polygon with respect to t° and ^ 'f'r+iix). 

For any negative rational number A, let P^ (—)(?/) S Vr[y] denote the residual polynomial in 
order r, with respect to t", 4>r{^)i and let R\{—){y) S Vr[y] denote the residual polynomial 
in order r + 1, with respect to t, (j)r+i{x), A. 

Proposition 3.5. Let t'^ G tr-i(/) be a non-complete type of order r — 1, and let t ~ 
{t^; —hr,y—rj) ^ tr{f) be a branch of t*^ such that er = fr = I ■ Let (j)r+i{x) be a representative 
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oft, and let 4>'j.{x) — (f>r-i-i{x) be the same polynomial, considered as a representative o/t*^. 
Let P{x) G 'Zp[x] be a nonzero polynomial. Then, with the previous notations, 

a) iN')'^^iP)^niN-+,{P)), 

b) mdt(P)=indt^,^,(P),. 

c) There exists a nonzero constant e e that depends only on t^, such that, for any 
A = —h/e, with h,e positive coprime integers, we have Rx{P){y) = €^R\_f^ {P){e'^y), 
where s is the initial abscissa of the X-component of N~j^^{P) [HNl §1.1]. 

Proof. We shall denote by Vr+i the p-adic valuation attached to t, and by Vr the p-adic 
valuation attached to t". Consider the (/i),.+i-adic development of P{x), which is simultaneously 
its 0j,-adic development: 

P[x) ^^ai{x)(l)r+i[xY ^^ai{x)(t)'^[xy. 

0<i 0<i 

For any < i, denote Ui = Vr+i{ai(j)\,j^i), u[ — Vr{ai{(j)'^Y), so that the points (i, Ui) determine 
the Newton polygon Nr+i{P), and the points determine the Newton polygon N^{P). 

Since degOi < to^ = Wr+i, 

Vr{ai) = ei • ■ ■ er-iv{ai{0)) = ei • • • e,._ierw(ai(6')) Vr+iioi), 

where 6 is any root of /t(a;) in Qp. By [HNl Prop.2.7+Thm.2.11], 

Vr+i{(t>r+l) — fr^rVr+li't'r) — Vr+li'Pr) = e^Url^r) + = '>^r{(t>r) + h^ — Wr(0'r) + K- 

This proves a), because Ui — = Wr(ai(0r)*) + "ihr = + ihr. 

Item b) is an immediate consequence, because Ti. transforms the horizontal line that passes 
through the last point of iV^j^(P) into the line L^hr of slope —hr that passes through the 
last point of N^{P) (cf. the figure below). 

Let us prove c). The definition of the residual coefficients and the residual polynomials is 
given in [HNl Defs.2.20-2.21]. Denote TV' = (iV')J?' (P). To every integer abscissa, < i < 
i{N'), one attaches a residual coefficient Ci of N^_^j^{P), and a residual coefficient c^ of N' , 
given by 

^. ^ f Zr''^'^^ Rr{at){zr), if (i , Mj) lics OU A^jT^^ (P) , 
' \ 0, otherwise. 

^ J zlZi'^^^ Rr-i{ai){zr-i), if lies on N' , 

' 1 0, otherwise. 

By a), the points {i,Ui), lying on the respective polygons have the same abscissas. 

Suppose that i is such an abscissa. For j — r,r — 1, denote by Sj{ai) the initial abscissa 
of the Aj-component of Nj{ai) [HN[ §1.1]. Since = 1, we can choose £r — 0. Since 
deg(aj) < rur, the polygon Nr{ai) is reduced to the point (0,Wr(aj))- This implies that 
tr(i) — (sr(ai) —£rUi)/er = 0] also, Rr{ai){y) is a constant, equal to z^J'Si'"^^ Rr-i{oi){zr-i). 
The exponents t'j,_i(i), and tr-i{0) are given by 

= (Sr-l(ai) - ir-lVr{ai{(f)'^y))/er-^l, ir-l(O) = (sr-i(aj) - er-iVriai))/er^i. 

Hence, c, = e'c^, where e = . Since Rx{P){y) = Cs + c,+,y + - ■ ■ + Cs+dey^ 

Rx-hAP)iy) = < + c'^+eV + ■■■ + <+de/, we get Rx{P)iy) = e^R'^_^^{P)ie^y). □ 
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Proposition 13. 5i applied to P{x) — f{x), shows that (p'^ix) is a better representative of 
t*^ than (j>r{x), in what the analysis of the branch t concerns. Actually, we have proved 
something stronger: the information obtained by applying to t the Basic algorithm in order 
r + 1, is exactly the same information obtained by applying the Basic algorithm to t° in order 
r, as long as we take (ji'ri^) ^ ^ representative, we analyze (^')r''(/) instead of the whole 
{N')~{f), and we replace indt(/) by ind^o"^^, (/). By "to obtain the same information" we 
mean to obtain the same number of new branches, the same (decreased) value of Wr+i(/) for 
each of them, and to cover the same distance to the end of the analysis of the branch t. 

Moreover, let A G be the slope of a side of N^_^^{f ), xj;{y) a monic irreducible factor of 
R\{f){y), and t" = (t°; —hr, ipr+i', A, ip{y)) the corresponding branch of t of order r + 1. By 
Proposition [SiSl this branch mirrors a branch of order r, t' — ((t°)'; A — hr, c'^°^^'ip{c~^y)), of 
t", with respect to the choice of (j>'j.{x) as a representative. 

Corollary 3.6. The types t' and t" are f -equivalent. 

Proof. If A = hr+i/cr+i, with /ir+i, Gr+i positive coprimc integers, then X — hr = (ft-r+i — 
Cr+ihr) / er+i, with coprime numerator and denominator; thus, e\.^i — e\ = e^+i- Also, 
fr+i = /r = degV', so that to*_|_2 = Br+i deg V^degc/iy+i = By c) of Proposition [3?5l 

applied to P{x) — ft'{x), we have w*+2(/t') = ^t+iiff), and 



This shows that ft'{x) has type t"; since deg/f — deg/f, we have necessarily /f ~ ft". □ 

This observation leads to an important optimization of the Basic algorithm. Whenever 
we apply the Main loop to a type and one of the outputs is a non-complete branch t = 
(to; -hr, y - 77), with e.^ = /,. = 1: 

(1) we replace t by the type t° itself, but taking (j)r{x) = 4>r+i{x) as a new representative, 

(2) we store the cutting slope hr as one of the data of the last level of t°, so that when 
the turn comes to apply the Main loop to the new t", only the sides of slope less than 
—hr of {N')~{f) will be analyzed. 

We call this a refinement step, and we use the term Montes ' algorithm to name the al- 
gorithm that is obtained from the Basic algorithm by applying a refinement step to every 
branch with e^/r — 1. By Corollarv l3.6[ every future branch t' of t is replaced by an equiva- 
lent branch t" of the new t°, so that the two algorithms are equivalent. By Proposition 
13. 5( the distance to the end of the algorithm covered by one application of the Main loop of 
Montes' algorithm is measured by ind|.o' (/). 



deg/t' = ml^^ujl^^ift') 



= m. 



t" t" 



iff)- 
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However, the refinement steps cause a strong diminution of the complexity. In fact, passing 
from order r to order r + 1 introduces a new level of recursivity in the basic tasks of the Main 
loop. Therefore, in Montes' algorithm the same information is obtained by working in a lower 
order. For instance, if in the Basic algorithm we find a branch of order r + n, with n successive 
levels with Cr+ifr+i = 1, for < i < n, 

t = (0i(x); Ai, 02(a;); ■ • ■ ; \r-i,4'rix);-hr, (j)r+l{x); • • • ; -/l^+n-l, 4>r+nix)). 

Starting with the trunk ((/)i(a;); Ai, (t>2{x); • ■ ■ ; Ar—i, (t>r{x)), we reached t by applying the Main 
loop n times in orders r,r + 1, . . . ,r + n. If we refine, this type collapses to 

t' = {(pi{x);\i,<j)2{x); - ■ ■ ;Xr-i,(j>'r{x)), 

with = (j)r+n{x), and we reach t' by applying the Main loop n times too, but always in 

order r. 

In order to homogenize the flow of the algorithm, we introduce a variable that stores 
an integer value at each r-th level of each type. Initially, it is given the value zero, and it is 
changed to — hr if we fall in a refinement step. This allows us to use a general Main loop, 
which presents only two differences with respect to the Main loop of the Basic algorithm: 

• in step 1 only the sides of slope less than —Hr of the Newton polygon of f{x) are 
analyzed, 

• in step 4, after computing the representative 0r+i of a non-complete new branch of 
order r, we proceed as follows: if e^/r > 1, we take the new branch as one of the 
output branches, with the value Hr+i = 0. If e^/r = 1, we take the input type as one 
of the output branches, with Hr = hr, and — (pr+i as a new representative. 

We can interpret also the refinement steps as a search for the optimal representatives. The 
search is performed by applying the Basic algorithm and not enlarging the types till a branch 
with Crfr > 1 is found. 

Summing up, Montes' algorithm has the same number of iterations than the basic algo- 
rithm, but a much lower complexity. It works only with optimal representatives, and it works 
always at the minimum order possible till a new optimal representative forces to pass to a 
higher order. 

In spite of this apparent strong optimality, one could speculate on an improvement based 
on a more intelligent way to pass from an optimal type of order r to an optimal type of order 
r + 1. The search for an optimal representative is done by blind lifting of certain polynomials 
over finite fields to Z, and a blind application of the Basic algorithm (without raising the 
order). This is extremely fast in practice, but there could exist a more direct way to obtain 
the next optimal representative. 

3.3. Computation of the index v^rith Montes' algorithm. Denote by t°P*(/) the set of 
types of order r that are produced by Montes' algorithm. The sets t°P*(/) are quite different 
from the sets tj.(/) produced by the Basic algorithm, which were crucial in the definition of 
indr(/) and the proof of the Theorem of the index. We need to compare in some sense these 
two types of sets. This is provided by Proposition 13.81 below, which is similar to Proposition 
I3.5( but going in the opposite direction. 

Let's go back to the situation of Corollarv l3.6l Suppose that the Basic algorithm is working 
with a type t° of order r — 1, and it finds a branch of order r + 1 (we denote it now by t 
instead of t"): 

t = {i° ; ~hr, <j>r+i; Xr+l,A+l{y)) ^ i'r+l{f), 

as the result of two consecutive enlargements of to, and we have e*/* = 1. Let c = e^^+'^ be 
the constant of c) of Proposition 13.51 and consider the branch of order r of t", 

t' = ((t")'; a;, V;(2/)), a; = A,+i - ^PUv) = c^'+^^.+i(c-iy). 
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By Corollarv l3.61 t' is /-equivalent to t. 

Remark 3.7. Any representative (j)'{x) of t' can be taken too as a representative of t. 

In fact, along the proof of CoroUarv 13.61 we saw that m*_|_2 — fn^+i, so that (j>'{x) has 
the right degree; thus, it is sufficient to check that (jJ*_,_2 (</"') = 1? s-nd this is given by c) of 
Proposition 13.51 applied to P{x) = 4>'{x). 

Let Nr+2{—), denote the principal Newton polygon of order r + 2, with respect to t and 
(l){x) :— (t)'{x). Let {N')r+i{—) denote the principal Newton polygon with respect to t' and 
(p'ix). Denote 

F, := F*; = F* = F*+i; F := F*+2 = ¥r[y]/A+iiy) = Vr[y]/^r+iic-'y) = 

For any A = ~h/e, with h, e coprime positive integers, let Rx{—){y) e V[y] denote the residual 
polynomial in order r + 2, with respect to t, (p{x), A, and let )(y) S F[y] denote the 

residual polynomial in order r + 1, with respect to t', (l)'{x), A. 

Let us write e^+i — e.^r+i = , hr+\ — For the type t' we have ^* /i* — (fj,)* e* = 1. 

Since h\ ~ ^r+i — ^r+ihr, this can be written as 

For the type t we have — (^^+i)*6r+i — 1- Therefore, we can choose = , 

Proposition 3.8. Let P{x) G Z[x] be a nonzero polynomial. Suppose we choose £*^i = £* . 
With the previous notations, 

a) Nr+2iP) = iN')r+i{P). 

b) indt(P) =indt'(P). 

c) There exists a constant r e F*, depending only on t' , such that for any A — —h/e, 
with h,e coprime positive integers, we have R\{P){y) — T^R'-)^{P){T~^y), where u is 
the ordinate of the initial point of the X-component of Nr+2{P) ■ 

Proof. For any polynomial a{x) G Z[x], 

^^+2(a)/e*+i = ■'^^+2(«)/er-+i, 'y*+i(a)/e* = i;*+i(a)/er+i, 

are the ordinates at the origin of the hues L\{Nr+i{a)), Lx-hr{{^')r{a.)), respectively |HN1 
Def.2.5]. These two ordinates at the origin coincide by a) of Proposition 13.51 so that w*^2 = 
This proves a), and b) is an immediate consequence. 
Consider the 0-adic development P{x) = X]o<i '^«(^)'?^(^)*' ^^'^ denote Ui = v^_^_l{ai(j)''') ~ 
vl_^_2{ai4>^) , for all i > 0. Let {ci}i>o be the residual coefficients of Nr+2{P), and {c^}z>o be 
the residual coefficients of {N')r+i{P). Since the two polygons are constructed from the same 
set of points {i,Ui) of the plane, we have q = if and only if c[ ~ 0. Let i be an integer 
abscissa such that the point {i,Ui) lies on Nr+2{P) ~ {^')r+i{P)), so that qc^ 0- f^is 
case, we have by definition. 

By definition, z^+i = y (mod ?/'r+i(2/)), and z'^ = y (mod Tpr+i{c^^y)) in Fr[j/]; thus, czr+i = 
z'^. By a) of Proposition [X5l applied to P{x) ~ ai{x\ we have s\j^][ai) — s\ {ai); since i^+i = 
£*. by hypothesis, and e*^^ = e* , we get tr+i{i) — t'j.{i). Therefore, by c) of Proposition [531 

4/c. = (z:,/^.+i)*''-(^)e-^''('^-) = e-^^'«-. 
If (s,m) is the initial point of the A-component of Nr+2{P) = i^')r+i{P), and i = s + je, we 
have Ui^u- jh, so that Rx{P){y) = T''R\{P){T-^y), for t = e^'' . □ 
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Therefore, we are able to deduce from the opthiial types constructed by Montes' algorithm, 
relevant information about the general types that would be constructed by the Basic algo- 
rithm. In particular, by an alternative and iterative application of Propositions 13.51 and [?751 
all values of indt(/), for all t e tr{f), can be captured along the flow of Montes' algorithm. 

Remark 3.9. If at the end of step 1 of the Main loop of Montes' algorithm, we accumulate to 
a global variable the value ind^o'^(/), the final output of this global variable is ind(/). 

In fact, if the input type t" is the result of an ordinary enlargement, then hr = and 
ind^o''(/) = iiidto(/); if t'^ is the result of a refinement step then, by Proposition l3.5i ind^o (/) 
is equal to the /-index of the type that the Basic algorithm would have produced if we had 
not refined. Proposition 13.81 guarantees that the future development of Montes's algorithm 
after a refinement step, mirrors the future development of the Basic algorithm. 



In this section we compute generators of the prime ideals lying above p in terms of the 
output of Montes' algorithm: a list T = {tp^, . . . ,tpj}, of /-complete types with optimal 
representatives, which are in 1 — 1 correspondence with the prime ideals pi,...,pg of K 
dividing pl^K- We write ef. A*!, etc. to indicate that a datum corresponds to the type 
tp. Recall that e{p/p) = • • • and f{p/p) — ■ ■ ■ fri can be read in the data of tp. We 
choose a root g Q of /(x), and denote by 9'^ G Qp the root of /t^ {x), image of 9 under the 
topological embedding K ^ Kp. 

If t S T has order zero, and (f>{x) is a representative of t, then the corresponding prime ideal 
is generated by (p, (f>{x)) by Kummer's criterion. If t S T has order one and its truncation to 
of order zero has indto (/) ~ 0, then the program computes generators of the corresponding 
prime ideal by using Dedekind's criterion. 

From now on, we fix a type t = tp, corresponding to a prime ideal p that did not fall in 
those special cases. We omit the superscript ( )P for the data of t. Let r be the order of t. 
We want to compute an integral element a — ap £ satisfying 



so that the ideal p is generated by p and a. 

Let us first construct an element P — f3p G K such that Wp(/3) = 1. To this end we compute 
first a representative (jjr+iix) of t. Since t is /-complete, ujr+i{f) = 1 and 4>r{x) ^ f{x). 
The Newton polygon 7Vr+i(/), with respect to t and (j)r+i{x), is one-sided of length one, and 
integer slope —H^ where H is the height of the side. By the theorem of the polygon, 



the last two equalities by [HNl Thm. 2. 11, Prop. 2. 7]. On the other hand, ujr+i{<l>r) = |HN| 
Prop.2.15], and [HNl Prop.2.9] shows that 



we have Vp{j3) = H. Thus, our aim is to find a kind of "worse possible" representative (jir+iix) 
of t; that is, one satisfying H = 1. To this end, we compute a blind (pr+iix). If = 1 we 
are done; if 7? > 1 we use a subroutine based on [HN[ Prop. 2. 10], to construct a polynomial 
P{x) € X[x\ with the following properties: 



4. Generators of the prime ideals 



wp(a) = l; Uq(a) = 0, Vq I p, q 7^ p; 



vi{a) > 0, V[tp, 



Vp{(t)r+l{9)) =e(p/p)z;(0r+i(6lP)) = Vr+l{(j}r+l) + H 

= erfrVr+l{4'r) + H ^ fr{e.rVr{4>r) + hr) + 




degP<mr+i, Vr+i{P) 



tV+l(<^r+l) + l, Rr{P){y) = l. 
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The point is that (pr+i '■= (j^r+iix) + P{x) is another representative of t, and it has H = 1. 
In fact, deg4>r+i = deg0r+i and Wr+i(0r+i) = t^r+i(0r+i) = 1, by |HN[ Prop. 2. 8], so that 
4>r+i is of type t and it 4>r+i is a representative of t. Now, since degP < rTi^+i, we have 

e(p/w e{p/p) e(p/p) 

Thus, v{^r+i{On) = <Pi0n) = K-+i(</'r+i) + l)/e{p/p), and 

Our next step is to compute the values Wq(/3), for aU other primes q lying above p, c{ =/= p. 

Definition 4.1. We say that tq dominates t, and we write tq > t, if tq is a branch of (t)r-i 
originated from a side o/A't,0^(/) of slope A < A^- In this case we denote Ap = A and we call 
it the dominating slope of tq over t = tp . 

Proposition 4.2. Let q 6e a prime ideal of K lying above p, q 7^ p. Let s be the order of tq. 
Then, 

' erfM---e^){\l-Xr). if tq >t, 

0, otherwise . 



Proof. Let tq be minimal with the property {tq)ra 7^ t,.Q. Necessarily rg < min{r, s}, because 
the types t, tq are complete. Let us deal first with the case ro < r. Since {tq)ra-i = tro-i, 
for some primitive choice 4>ra{x) of a representative of this type, the Main loop produced 
(at least) two different branches, that later developed to produce the types t, tq. Some 
of these branches might have been refined, causing a change of this representative at level 
ro — 1; however, arguing with these primitive branches if it were necessary, we can assume that 
(prg — 4>f.g = (fira ■ This might change the values of the data of the r-th level, but this is not 
relevant in our arguments. After our assumption, vf = v'^, Nf{—) = N^{—), for 1 < « < ''o- 
We claim that 

(4.1) v^icl}r+i) = vU(b''/^''), foralll<z<ro + l. 

Let us show this by induction; clearly, Vi{(t)r+i) = = w5'((/)^'^''^''), because both polynomials 
are monic. Suppose that (|4.ip holds for some 1 < i < rp, and let us show that it holds for 
i + 1. By the definition of v'^_^^^, we need only to show that 

7V^(0.+i)=iV^(0^'-/'-). 

Since i < ro, it is sufficient to check that Nf{(t)r+i) = Nf {(p'^''^''). Now, these polygons are 
both one-sided of slope A^, and have the same length because the two polynomials, (j)r+i, 4''?'^"' i 
have the same degree. Finally, the two polygons have the same end point by the equality of 
(|4?T|) for our z, and by |HN[ Lem.2.17]. This ends the proof of (|4T|) . 
We claim now that 

In fact, the polygon N^^{(f)r) = Nf^{(l)r) is one-sided of slope \r„. If X^^ ^ A,.o, then R^^[(f)r) 
is a constant and ujl^^j_-^^{<j)r) = 0. If — A^o, but i/'ro V'roj then i?J!^(0r) — Rra{4>r) is a 
power of ipro up to a multiplicative constant, and t/j^^ \ R'^^{(j)r), so that uj^^^j^{(j)r) = too. 
The same argument works for (pr+i- 

Finally, (gT]) and gS]) show that v{(j)r+i{e'^)) = v{(j)r{e'^Y^-f-). Therefore, Vq{(j)r+i{e)) = 
Vq{cj,r{OY^-f^-), and Vq{l3)^Q. 

Assume now ro = r. The Main loop applied to {tq)r-i = t^-i and the representative 
<j)r — (t>r produced a complete branch t, and (at least) another branch of order r 

t° = iM^);--- ;A,_i,(/.,(a:);A°,V°(y)), 
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whose further development produced the type tq. Let us denote by v^^i, ^r+ii the data 
attached to this type. Since O'^ is the root of a p-adic polynomial of type t[|, the Theorem of 
the polygon shows that, 

(4.3) viMO'')) - iVri<f>r) + |Aj!|)/ei • • • er-l. 

Let us compute now v{(j)r+i{0'^))- The Newton polygon N^{(j)r+i) = N^{4ir+i) has length 
Cr/r, slope Ar, and the ordinate of the last point is Vr{4)r+i) = Vr{4'P'^'')- Consider a line of 
slope Aj! far beyond the polygon, and move it upwards till it touches it; let L be this line of 
slope A" that first touches the polygon, and let H be the ordinate of L at the origin. We 
distinguish two cases, according to A° < A^ or A° > A^. In the first case, A° = Ap is the 
dominating slope of q over p. 





A^ <^ \r A^ ^ A^ 

Arguing as in the case tq < r, we get = in both cases, and 

(4.4) v{<t>r+i{B'^)) = f°+i(0.+i)/ei • • • e,_ie° = H/e^ ■ ■ ■ e,_i, 

the last equality by the definition of f^^^. The figures show that H — Vr{(j>r''-^'') 

H' = 



H\ with 



G-f fr I Xf I , 

erfr\X% 



if AO < A,., 
if AO > A,.. 



Now, (jiral) and (|4r4|) . show that 



ei • • • Cr^l 



ei • • • er-i 



Therefore, Wq(/3) = if Aj! > Ar, and otherwise. 



ex-- ■ Cr-l 



-e./.(|A.|-|Ag|). 



□ 

For the maximal types with respect to the ordering ">", we take dtp :— f3p. For the rest of 
the types we compute recurrently: 



n 



These elements are not far from generating the pi, since: 

I if q = p, 
otherwise. 



Unfortunately, they could be non- integral at primes of not dividing p Zk ■ This can be 
easily arranged; we write each ctp in the form cup = G{9)/b, with G{x) G Z[x] and & G Z 
coprime with the content of G{x): and we conveniently modify dtp into: 



an 



G(6i)/p''(''). 
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5. Computational issues 

Recall that Montes' algorithm is the optimization of the Basic algorithm of section [2T2] that 
results from the application of the refinement process of section [31 In this section we give a 
more detailed description of this algorithm and we discuss some computational aspects. We 
also include, as an extension of the algorithm, the computation of generators for the prime 
ideals as indicated in the last section. 

5.1. Outline of the algorithm. The primary goal of Montes' algorithm is the computation 
of ind(/) and the construction of a set T of /-complete types, that faithfully represents f{x). 

By the recursive nature of its construction, many of the types generated by the algorithm 
will share many of its levels, so that most of the computations necessary to enlarge them will 
be the same. Hence it is very convenient to organize their computation in such a way that we 
can take profit of as much previous computations as possible. The simplest way to organize 
the computation of types is to store all the types being built by the algorithm in a list, which 
we call STACK. Once a type is completed, it is moved from the list STACK to a second list called 
COMPLETETYPES, which when the algorithm ends contains the final list of complete types. 

The variable STACK, as its name suggests, is a LIFO stack, which in practice determines the 
flow of the algorithm: the main loop of the algorithm extracts the last type from the STACK 
and works it out to decide whether it is complete (in which case it is moved to COMPLETETYPES) 
or it originates a number of enlarged types, that will be added to the top of the STACK. The 
program finishes when the STACK is empty. 

We could think of the types being built along the algorithm as the branches of a tree. 
The root of the tree is a node corresponding to the input polynomial f{x). Every division of 
the branch into new subbranches is generated in every order by multiple sides of a Newton 
polygon of /(x) and by multiple irreducible factors of the residual polynomial of each side. 
The algorithm builds this tree of types from the topmost branch to the lowest one in every 
order. This strategy confers a certain ordering to the list COMPLETETYPES, which is useful 
later on for the computation of generators of the ideals. 

The computation of the p-index is performed along the construction of types: every time 
we analyze a Newton polygon with respect to a type t, we add the number ind^' (/), given in 
(f3?2|l . to the variable TOTALINDEX, whose final output is the value of ind(/). 

Once the algorithm has emptied the STACK, the algorithm is almost finished: it remains 
only to gather the information of every complete type to list the ramification indices and 
residual degrees of the prime ideals dividing jiIk- 

We now give a detailed outhne of Montes' algorithm, using standard pseudo-code. 

MONTES' ALGORITHM 

INPUT: 

- A monic irreducible polynomial /(x) G 7L\t\. 

- A prime number p G Z. 

OUTPUT: 

- The p-valuation of the index of ]{x) in TLk- 

- A list {(ei, /i), . . . , (eg, /g)} of pairs of integers describing the factorization of p1,K'- 

- A list of integral elements ai, . . . , Ofg G Im such that pi — pZ,K + chiLk 
INITIALIZATION STEPS 

1 Factor /(y) = "01(2/)°^ ' ' ' V's(2/)"° modulo p, with ^/'i(y) G Fp[j/] pairwise different monic 
irreducible polynomials. 
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2 Take monic polynomials 4>i{x), . . . , (j^gix) G Z[a;] such that 4'i{y) (mod p) — ipi{y). Compute 
the polynomial M{x) = (/(x) - • • • (l)s{xY=)lp. 

3 Initiahze empty lists STACK and CDMPLETETYPES, and set TDTALINDEX^ 0. 

4 FOR every polynomial (t)i{x) do 

(Dedekind's criterion) If ai = 1 or i!i{y) \ M{y) (mod p), output the ideal p — (p, (t)i{9)), 
with ramification index e{p/p) = ai and residual degree f{p/p) ~ deg Otherwise, add 
to STACK the extension t = {(j)i{x)) of the type of order zero determined by ilJi{y), and set 
Hi =0, as data of level one. 

MAIN LOOP 

WHILE the STACK is non-empty do: 

5 Extract the last type t° from STACK. Let ?■ — 1 be its order. 

6 Compute the Newton polygon N^''{f), formed by the sides of slope smaller than —Hr- 

7 Compute ind^" (/) using the formula ([3?2l) and add this number to TOTALINDEX. 

FOR every side S of N^-{f) do 

8 Set \f ^ slope of S. 

9 Compute the r-th order residual polynomial Rr{f){y) £ ^r[y]- 

10 FOR every irreducible factor ^'(y) of Rr{f){y) do 

11 Make a copy t of the type t°, and set ^'^(z/) ^ "0(2/) ■ 

12 Compute a representative (j)r+i{x) G oft. 

13 If ovA^{Rr(f)) = 1 (the type is complete), add (t; A*, ■(/'*(y)) to CDMPLETETYPES, and 
continue to the next factor of Rr{f){y). 

14 If degV' = 1 and A* g Z (the type must be refined), set (l)l{x) ^ (f>r+i{x), Hf. ^ A*, 
add t to the top of the STACK and continue to the next factor of Rr{f){y). 

15 (Build a higher order type) Add (t; A*, V'*(j/)) to the top of the STACK. 

END OF MAIN LOOP 
OUTPUT 

16 Print the p- valuation of the index of / in , given by the value of TOTALINDEX. 

17 For every type t in the list CDMPLETETYPES, output the ramification index and residual 
degree of the corresponding ideal p, given by: 

e(p/p) = e*...e*, /(p/p) = m^i* ' ' ' 
where m* — deg 4>\ , and r is the order of t 

EXTENSION: COMPUTATION OF GENERATORS 

All notations are taken from section [4] 

18 FOR every type tp in CDMPLETETYPES compute the element /3p = 4)r+i{9)/<j)riOy'-^'' ■ 

19 FOR every type tp in CDMPLETETYPES compute the element ap := /3p nt,>t "q^"''*^'^ 
where Vq{/3p) is given in Proposition |4?2] 

20 FOR every type tp in CDMPLETETYPES compute the element ap = G{x)/p"'^^\ 

To compute ap in step 19, it is necessary to slightly modify the algorithm in order to store 
all dominating slopes Ap. 
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5.2. Some examples. 

Example 1. Let us consider the irreducible polynomial 

f{x) :=x^^ - 588a;^" + A76x^ + 130095x^ - 172872a;^ - 12522636x^ + 24745392a:^ 

+ 486721116a;^ - 1583408736a;^ - 641009376a;2 + 10978063488a; + 59914669248, 

whose discriminant is 

disc(F) = 2^^ ■ 3^^ ■ t"^ ■ 79'' • 14159^ • 644173^ • 3352073^ . 

We apply the algorithm to find the decomposition of the prime p = 2 in the ring of integers 
Zif of the number field K — Q{9) generated by any root of the polynomial f{x). Since 

f{y) = iy + l fy' (mod 2), 

we find two types ti, t2 of order zero, extended respectively to (j>\ (x) = a; + 1, = ^- The 

Newton polygon Ni^{f) has two sides, with slopes —3/2 and —1/2 respectively, which single 
out two prime ideals pi,p2, with e(pi/2) = e(p2/2) = 2 and /(pi/2) = /(P2/2) = 1. The type 
tpj dominates tp^ with dominating slope Xpl = —3/2. 

The Newton polygon A^*^(/) has again two sides, with slopes —1 and —1/2, and residual 
polynomials Ri.i{f) = {y + 1)'', Ri,2{f) = {y + 1)^, respectively. Hence, the type t2 yields 
two types t2,i,t2,2 of order one. The first one must be refined and the second one must be 
enlarged to an order 2 type. To refine t2,i, we take the new polynomial (f>i'^{x) ~ x + 2. 
The corresponding Newton polygon has only one side with slope smaller than —1; the slope 
is —3/2 and the residual polynomial {y + 1)^, so that this type must be enlarged too, to an 
order 2 type. After computing their respective representatives, we have now two extended 
types of order one, ready to be enlarged to order 2: 

t2,i = (a; + 2;-3/2,(a; + 2)2 + 8), i?i = -1, 
t2,2 = (.x;-l/2,a;2 + 2), Hi = 0. 

The Newton polygon N2^'^ (/) has a unique side with slope —4 and residual polynomial (y+l)^, 
so that this type can be refined. We take, for instance, (p^'^ {x) ^ {x + 2)^ + 40: 

t2,i = (.T + 2;-3/2, (a; + 2)2 + 40), i?i = -1, H2 = -4. 

The new polygon iV2^'^(/) has two sides with slopes —9 and —5, originating two new prime 
ideals P3,p4 dividing 2Zk, with e(p3/2) = e(p4/2) = 2, /(P3/2) - e(p4/2) = 1. The type 
dominates tp^ with dominating slope Ap^ — —9. 

The Newton polygon N2^'^{f) has a unique side with slope —4 and residual polynomial 
{y + 1)^; thus, we must refine. Take 02^'^ (a;) — x"^ + 10, H2 = —4. The next Newton polygon 
has again a unique side with slope —5 and residual polynomial {y + 1)^. We refine again, 
taking (j>l''-^{x) = a;^ + 10 + 8a:. We get: 

t2,2 = [x; -1/2, x^ + 8x+ 10), Hi =0, H2 = -5. 

Now, N2^'^{f) has two sides with slopes —8 and —7, both with residual polynomial y + 1, 
thus giving two prime ideals p5,p6, with e(p5/2) ~ e{pe/2) ~ 2, /(P5/2) ~ e(p6/2) = 1. The 
type tp5 dominates tpg with dominating slope Ap^ — —8. 

Summing up, 2Zk = (pi • ■ -pe)^- The 2-index of /(x) isind2(/) = 33, and u (disc (i4')) = 18. 

Generators for the ideals pi can be determined by using the formulas of section|4l We find: 

. _ (g + 1)^ + 4(6) + 1) + 8 , _ (6> + 1)2 + 2(61 + 1) + 2 
"1- (9+1)^ ' (0 + 1)2 

(61 + 2)2 + 64(61 + 2) + 40 , _ (6* + 2)^ + 16(6* + 2) + 40 .4 



(6* + 2)2 + 40 ' (61 + 2)2+40 

6i2 + 4061 + 42 _ _ 612 + 2461+ 10 . 

' + ' " 6^ + 89 + 10 ' 
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ai = (261" + 6li" + 26i9 - 6^ + 2)/2; 

aa = (61" + 6li" + 6*9 - 36*8 + 46*^ + + 4) /4; 

as (13161" - 4746110 + 4486*9 + 526i8 + 30961^ - 3666*6 _ 2i66»5 + 2566l4 - 3646l3 

+ 1366*2 - 4966* + 32)/512; 
UA = (-276*" - 26^° ~ 806*9 ^^2618 + 196*^ + 106**^ + 726*5 _ 646*" + 766*^ 

-1046*2 + 486* + 32)/128; 
as = (-456*" - 106*10 - 186*9 + 646*^ - 76*^ + 506*0 _ 39515 _ qqqA ^ 325/ _^ 64)/64; 
ae = (336*" + 86*10 + 426*0 _ 4518 _ 53517 _ _ ^^q5 _^ ^^qa ^ 325/3 _^ 54512 _ 325i)/64. 

Example 2. Take p — 2 and consider the irreducible polynomial 

!{x) {x^+x + 5)50 ^ 289(2;3 + a; + 5)^5 + 2^^^ . 

The algorithm takes initially 4'i{x) = x'^ + x + 1, and finds a unique side with slope -2 and 
residual polynomial (y+l)50. A refinement leads to 01 (x) = x^ + x + ^>^ and a Newton polygon 
with one side, with slope —89/25 and irreducible residual polynomial + y + 1. Hence, in 
the number field K defined by any root of f{x), we have 2'Lk = P^^, with /(p/2) = 6. The 
2-index of the polynomial is 13011. While this computation is almost instantaneous, the 
determination with Pari of a 2-integral basis of K takes about 190 seconds, and needs an 
amount of 244 Mb of memory. 

5.3. Some remarks on the complexity. We have not developed a detailed analysis of the 
complexity of the algorithm, but the experimental results of the next section indicate that its 
running time is excellent. We now provide some arguments to explain this good behaviour. 

We saw in section [^751 that the number of iteration of the Main loop is bounded by ind(/) 
and that each iteration covers indt(/) steps from the total value of ind(/). One is tempted 
to conclude that the running time of the algorithm is linear on the discriminant of f{x), 
but this is not so evident, since the treatment of higher order types is much expensive than 
the treatment of low order types. However, this is balanced by two facts: on one hand, 
indt(/) is generally much bigger than one in each iteration; on the other hand, the higher 
the order, the smaller is uJr+l{f)^ and this invariant tells the number of coefficients of the 
<j)r-a.dic development of f{x) that must be computed, the length of the Newton polygon to be 
analyzed, and it is an upper bound for the degrees of the residual polynomials. 

At least, it seems that for polynomials whose types have bounded order, the running time 
will be at most linear on the discriminant. In practice, the average running time is much 
smaller. 

On the other hand, the degree of the polynomials 4>\{x) appearing in an /-complete type t 
is a divisor of the product e(pt/p)/(pt/p)- Every time the type is enlarged, the degree of the 
last 4>k{x) is multiplied by the corresponding product e*/*. Hence, for a polynomial f{x) to 
have attached a type of a very high order, its degree must be really huge. This explains why 
the algorithm works well for polynomials of high degree: the maximum of the orders of the 
types of a polynomial grows slowly in comparison with the degree. 

The low memory requirement of the algorithm is another of its strong advantages: it is 
only necessary to store the polynomials (pkix) and tpk{x) of the types being calculated (and 
some Taylor expansions to gain efficiency) . This makes possible the treatment of polynomials 
of very high degree with scarce computational resources. 

The complexity of the computation of the generators is dominated by the inversion of (pr (0) 
in K, which is a hard task if the degree of 4>r{x) is large. 

5.4. Implementation of the algorithm. The first implementation of Montes' algorithm 
was programmed by J. Guardia in 1997, as a part of his Ph.D. It was written for Mathe- 
matica 3.0, and it included a specific package to work with finite fields, since that version of 
Mathematica did not carry such a package. It is still available on request to the author. Ten 
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years later, we started a collaboration to make a full upgrade of the algorithm, with many op- 
timizations both theoretical and computational, including a completely new implementation 
in Magma. 

The computation of generators for the prime ideals becomes a heavy time-consuming task 
if the p-adic factors of f{x) have large degrees. For this reason, our implementation skips 
this calculation by defect. If the user wants to compute the generators, a Boolean variable 
GENERATORS has to be given the value "true" . 

As memory requirement is not a constraint for the implementation, the program stores 
some intermediate results (mainly (/)-adic expansions) to gain speed. The main data type 
used by the program is a specifically designed record which contain all the relevant data of a 
type in a given order. To avoid massive replication of the types being computed, must of the 
routines access them by memory address. 

The program also includes a number of routines to construct types and polynomials with a 
prescribed set of types. The program and its documentation, which includes all the examples 
presented in this paper, can be downloaded from the web page 

www-ma4 . upc . edu/ ^guardia/MontesAlgorithm. html. 

Any comment on the program will be welcome. 



6. Some heuristics on the complexity 

We dedicate this section to illustrate the performance of (our implementation of) Montes' 
algorithm with several polynomials chosen to force its capabilities at maximum in three di- 
rections: polynomials with a unique associate type of large order, polynomials which require 
a lot of refinements, and polynomials with many different types. We have also included some 
test polynomials found in the literature. 

All the tests have been done in a personal computer, with an Intel Core Duo processor, 
running at 2.2 Mhz, with 3Gb of RAM memory. The reader willing to check these results on 
his/her own can obtain the Magma code to generate these polynomials from the web page of 
the program. 

Example 1: Take p= 2. Consider the irreducible polynomials 

(l)l=X^ + + 2^; 
02 = + 24x<^i + 212; 

<f>3 = <Pi + 2^3 (x + 22)</,2 + 2^^xcPv, 

<t>5 = <i + 2^^4>Xh4>l + 2'^%a;(<^i +2^){4>l + 2^^<k2) + 2''h)<t>3 + ^'^x^i^Pl + 2^3)); 
cj>e = <Pl + 2^^Ms4 + 2^^^^((a; + A)cf>, + 2^)cj>lcf>r, 

<j>T = cj)l + 22351((^,03 ^ 2^^x{^, + 2^)^2)(f>4 + 2'''>^ {x<f>i(f>l + 225((x + 2^)<i>i + 2^x)<l>2M 
+23234(((^<^^03 + 225(x + 2^){^, + 2^)4>2)h + 2^°(x</.i + 2^''))<t>i 

+2i68((a; + 22)<^i+26a;)<^2)^^. 

08 = 275i5((((((a; + 4)01 + 26x)0i + 231x)03 + 239.T0103 + 2^0x02)01 

+2l04(((^<^^ + 28)0i + 225(x0i + 2\x + 24)))03 + 238((x + 4)(0i + 26)03 + 2^701 02))04 

+22O8(((x0i + 28)0i + 225x01 + 23i(x + 4))03 + 2^Vl<^2 + 2*^^(0:01 + 28)02))0i 

+2924(((((:c + 4)01 + 26x)03 + 225((x + 4)01 + 28)02)03 
+2i34((^, + 24(x + 4))0i + 225(01 + 24x)))0l 

+2l04(((a; + 4)0103 + 225x(0i + 26)02)03 + 264((x + 4)(0i + 2^)02 + 2270i))04 
+221O(((0, + 24x)03 + 223(x + 4)(0i + 2^)02)03 + 2^^(0102 + 2250i)))02)0g03 

+2206i8(((^^^<^3 + 23i(x + 26)02)03 + 2^%ct>^4>l + 2^'{<^, + 2^)))4>l 

+2104(((X + 4)0103 + 225x0102)03 + 270((x + 4)02 + 22l0i))04 

+2208((((x + 4)01 + 26x)03 + 225((x + 4)0i + 28)02)03 + 2™(x02 + 2l9(x0i + 28))))0i 
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+2i"4((((a; + 4)0^ + 26a;)02 + 233)03 + 264x(/.i<^2)<^4 
+ 22"8(a; + 4)(0i + 26)02^3 + 2249(01 + 2^x)<j)^)<j)l 

For each j, the corresponding polynomial (pj has a unique associate complete type of order 
j, so that in the corresponding number field Kj the ideal 2Zxj is the power of a unique prime 
ideal pj . The following table contains the degree and 2-indcx of the (f>k , the ramification index 
Cj and residual degree fj of pj and the time ti used by the program to compute them. The 
last column indicates the time t2 to run the extended version of the algorithm, which includes 
the computation of generators for the ideals pj. All the times arc expressed in seconds. The 
computation of the generators in the last two rows was stopped after 24 hours of running 
time. 



(j)j 


deg 0j 


ind(0j) 




/. 


tl 


t2 


<t>l 


2 


2 


1 


2 


0.00 


0.00 


4>2 


4 


16 


1 


4 


0.01 


0.01 




16 


360 


2 


8 


0.01 


0.01 


04 


32 


1544 


2 


16 


0.01 


0.016 


4>-> 


96 


14616 


2 


48 


0.08 


0.6 




576 


537120 


6 


96 


0.70 


406 


07 


1152 


2153376 


12 


96 


4.0 






6912 


77673504 


36 


192 


787 





Example 2: Let f^{x) — {x'^ +a; + 1)^ — |j2fe+i^ with p = 1 (mod 7) a prime number. When 
we apply Montes'algorithm to factor the ideal pl^K, we obtain two types of order zero with 
liftings 01 (a;) e Z[x] of degree one. For both of them the Newton polygon has only one side, 
with slope —1 and end points (2, 0) and (0, 2), and the residual polynomial is the square of a 
linear factor. After approximately 2k total refinements, both types become /'^-complete. The 
ideal pl^K splits as the product of two prime ideals with ramification index 2 and residual 
degree 1, and the p- index of f*'{x) is 2k. 

This is almost the illest-conditioncd quartic polynomial for the algorithm, since the index 
of every type is increased a unit per refinement in general, and the total p-index of f^{x) 
is 2k. Thus, the program has to make about 2k iterations of the main loop. Numerical 
experimentation shows that even in this worst case the running time of the algorithm is very 
low. In the following table we show the running time of the programm for different values 
of k and p. As before, ti is the time in seconds to compute the index, residual degrees and 
ramification indices, and t2 is the time to compute also the generators for the prime ideals. 



p 


ind(/'') 


tl 


t2 


P 


ind(/'=) 


tl 


t2 


7 


1000 


0.57 


0.62 


43 


10000 


229 


237 


7 


2000 


1.95 


2.1 


103 


10000 


324 


334 


7 


4000 


8.7 


9.2 


1009 


1000 


2.1 


2.6 


7 


8000 


44.7 


46.2 


1009 


2000 


10.8 


12.3 


7 


16000 


245 


250 


1009 


4000 


58 


62 


7 


20000 


436 


444 


10^^ + 9 


1000 


10 


12.7 


13 


1000 


0.75 


0.85 


10^ + 9 


2000 


57 


66.5 


13 


2000 


2.9 


3.1 


10^ + 9 


4000 


313 


341 


13 


10000 


131 


135 


10«9 + 9 


100 


2.8 


4.8 


19 


10000 


158 


162 


10«9 + 9 


200 


8 


13.5 


31 


10000 


198 


205 


10^9 + 9 


400 


29.4 


48 


37 


10000 


214 


221 


lO'^y + 9 


1000 


221 


308 
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Example 3: Take p = 13. We now consider a polynomial with several different types. Let 

(j)i{x) = x2 + 132x+ IS'' • 3; 
</'2(x)=0i(x)3 + ((13i«.2)); 

<t>3{x) = <^2(x)io + 13S9(a; + 132)02(x)5 + 13™^i(a;); 

(^4(x) = 03(x)2 + 13248(12(a; + 132)0i(a;) + 138)02(x)6 + 13^35 • 120i(x)202(x); 
i 

- Y[ + k) + 13^"°", j - 0, . . . , 12; 

fe=0 

In the number field Kj defined by fj{x), we have the factorization 

li'LL=p\---p]. /(p,/13)-24. 

Each prime ideal comes from a different order 4 type. The 13-index of fj{x) is 21576j. The 
times to compute this index and the factorization of 13 in the fields Kj are shown in the table 
below. 



j 




ind(/,) 


ti 


1 


120 


21576 


0.08 


2 


240 


43152 


0.3 


3 


360 


64728 


0.9 


4 


480 


86304 


2.3 


5 


600 


107880 


4.4 


6 


720 


129456 


8.2 


7 


840 


151032 


13.3 


8 


960 


172608 


20.3 


9 


1080 


194184 


29.4 


10 


1200 


215760 


40.6 


11 


1320 


237336 


55.2 


12 


1440 


258912 


72 


13 


1560 


280488 


92 



The computation of the generators for the number field defined by polynomial fi{x) took 
20 seconds, for the number field defined by f2{x) lasted 61/2 hours. The computation for /a 
exhausted Magma's virtual memory due to a coefficient explosion in the computation of an 
extended gcd. 

In this example one cannot expect a linear behaviour of the time versus the number of types, 
because the addition of more factors to the product defining the fj{x) implies a significant 
growing in the size of the coefficients of the polynomial, which has a certain impact in the 
running time of the algorithm. 

Example 4: We applied the algorithm to the list of 32 polynomials /i, . . . , f^2 appearing 
in [FPRO2I, appendix D]. The total running time for altogether was less than 0.2 seconds. 
We then applied the algorithm to the polynomials Fi ~ ff +pl°'^^, where pi is the prime 
specified in loc.cit. for every polynomial. In the table below we display the index of these 
polynomials and the running times of the algorithm. As before, ti denotes the time in seconds 
to determine the index, the residual degrees and the ramification indices, and t2 is the time 
of the extended algorithm that includes the computation of the generators of the ideals. 
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/ 


P 


ind(/) 


ti 


t2 


/ 


P 


ind(/) 


ti 


t2 


Fi 


2 


2502150 


2.45 


3.1 


Fn 


2 


1571054 


2.9 


7 


F2 


2 


1481141 


2 


2.7 


Fis 


7 


7331055 


77 


93.9 


F3 


3 


2570992 


6.6 


8.2 


Fig 


71 


187219 


116.1 


249.3 


Fi 


3 


1177569 


4.6 


7.5 


F20 


3 


287752 


15.4 


23.2 


F5 


2 


2502505 


1 


1.5 


F21 


5 


10117231 


73.8 


81.2 


Fe 


2 


2137558 


1.6 


2.4 


F22 


3 


5194476 


18.4 


23.6 


F7 


2 


2751159 


3 


4.7 


F23 


3 


2888852 


15.7 


23.6 


Fs 


5 


1646099 


13.81 


20.5 


F24 


2 


2901708 


6.2 


9.5 


F9 


2 


1672713 


2 


3 


F25 


47 


2612660 


253.5 


636 


Fio 


1289 


1500768 


117 


234 


F2& 


61 


4257732 


158 


192 


Fn 


2 


2629928 


3 


4.1 


F27 


2 


5925350 


7.7 


9.4 


F12 


3 


5895414 


20 


22.6 


F2S 


3 


5720164 


5 


7.1 


F13 


11 


1810788 


35 


64.6 


F2g 


3 


7826660 


15.8 


23 


Fu 


17 


1618581 


31.7 


58 


-F30 


2 


39363539 


14.7 


26.3 


Fl5 


2 


7744913 


4.9 


6.1 


F31 


2 


40933692 


62.4 


73 


Fig 


2 


3808558 


4.3 


6.9 


F32 


2 


17097775 


82.7 


132 
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